Steady-State Properties of a Totally Asymmetric Exclusion Process with Periodic 

Structure 
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We formulate and analyze the steady-state behavior of totally asymmetric simple exclusion pro- 
cesses (TASEPs) that contain periodically varying movement rates. In our models, particles at a 
majority sites hop to the right with rate pi while particles occupying a periodically arranged set of 
sites move to the right at rate P2- A number of approximate mean field approaches are used to study 
the steady-state currents and bulk densities of this model. While exact solutions are not found, the 
mean field approaches provide results that show good agreement with data derived from extensive 
Monte-Carlo simulations. 
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I. INTRODUCTION 



As a canonical model of one-dimensional transport, the 
Totally Asymmetric Simple Exclusion Process (TASEP) 
is a topic of current interest. One of the few solvable 
non-equilibrium models, the TASEP and its solutions 
have been extensively investigated 0, 0, El El 13] as a 
model for numerous one-dimensional transport processes 
including ribosome motion H, IE j P ore transport @ , 
and traffic flow @, E3- The traditional TASEP model 
consists of a finite lattice with open boundaries. Parti- 
cles are inserted into an empty site at the leftmost end 
of the lattice at a rate a, and once in the lattice the 
particles move to the right at a rate p. Each motion 
of a particle within the lattice moves that particle ex- 
actly one lattice site to the right, and particles can move 
only if the site immediately to the right is not occupied 
by another particle. Upon reaching the rightmost lat- 
tice site, particles are removed from the lattice at a rate 
/3. Exact solutions for the steady-state particle currents 
(J), and densities (a), are available through mean field 
approaches 0,0, matrix product methods Q, and recur- 
rence relations The latter two approaches solve the 
model completely, providing exact expressions for density 
correlation functions of arbitrary order. The TASEP so- 
lution itself yields a phase diagram with three phases. At 
large values of the injection (a) and extraction (/3) rates 
the system is dominated by the rate of particle hopping 
and is in a maximal current phase. At small values of a 
and /?, the system is found in a low-density entry limited 
phase, and a high-density exit-limited phase respectively. 
Within each of these two boundary-limited phases lie a 
pair of sub-phases each with a distinct density profile. 
The fundamental form of this phase diagram has proven 



to be very robust and extending the TASEP to include 
particles that occlude more than one lattice site or 
backwards particle motions , has not altered the phase 
diagram significantly. Nonetheless, these extensions have 
facilitated the use of the TASEP in modeling a wide va- 
riety of physical processes. In many phenomena of inter- 
est however, the assumption of a single internal hopping 
rate p implicit in the normal TASEP does not adequately 
capture the full character of the transport behavior. One 
system that has been the subject of intense interest is 
the kinesin motor protein and related molecules ^jj . Ki- 
nesin, and kinesin-like molecules, form an integral part 
of the cellular transport system, shuttling protein-filled 
lipid vesicles to various parts of the cell. Traveling along 
tubulin polymers, the motion of the kinesin molecule is 
driven by a complex conformational change producing 
steps approximately 8nm in length |12|. The exact step- 
ping mechanism of the kinesin family of motor proteins 
is currently the subject of speculation, however one of 
the lea ding candidate models is the so-called " inchworm" 
model 13]- The inchworm model proposes that a sin- 
gle 8nm step is composed of two sub-steps; in the first 
sub-step one kinesin head group moves from its current 
position to the next available binding site on the tubu- 
lin filament. This sub-step is immediately followed by a 
dragging motion where the lagging head group is pulled 
up to essentially the same position as the lead group. 



a 



f 1 

1- k 






\ M 


■□□1 


• 


• 


m 


won 


W 




4 


p 2 




Pi 



'Electronic address: glakatos@physics.ubc.ca 
t Electronic address: tomchouOucla. 
* Electronic address: tolya@rice.edu 



FIG. 1: The two-rate totally asymmetric exclusion model. At 
every T lattice sites, the particle movement rate is p2. At all 
intervening sites, the movement rate is p\. All other aspects 
of the model are identical to those of the standard TASEP. 

While TASEP models have been successfully applied to 
model molecular motor dynamics in the past |14j , apply- 
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ing the TASEP to the two-part motion of the inchworm 
model requires that the TASEP be extended. As the two 
sub-steps in the inchworm model are not expected to oc- 
cur at the same rate we can expand the applicability of 
the TASEP to these types of systems by introducing two 
internal hopping rates pi, and P2 (cf. Fig. ^| into the 
TASEP model. We assign these hopping rates to the 
lattice sites in an alternating fashion where odd indexed 
lattice sites move particles to the right with rate p\ , while 
even indexed sites move particles with rate p 2 - Using this 
arrangement a single (pi,p 2 ) unit would comprise, for ex- 
ample, a 8nm kinesin step with the p\ step representing 
the movement of the lead head group, and the p 2 step the 
movement of the trailing group. Generalizing this two- 
rate TASEP model, we also consider the case where the 
Pi sites are separated by an arbitrary number (fc) of p\ 
sites. Applying mean field methods we develop accurate 
approximations to the steady-state current and density of 
the dual-rate TASEP in the maximal current phase, en- 
try limited, and exit limited regimes. We also conjecture 
that, as was true with many previous extensions to the 
standard TASEP, the dual-rate extension does not funda- 
mentally alter the form the TASEP phase diagram. The 
accuracy of our mean field methods are verified through 
extensive Monte Carlo simulations. 



II. MEAN FIELD THEORIES 

A. Simple Mean Field Methods 

We begin by considering the dual-rate TASEP in the 
limit of large a and 0. As is the case in the standard 
TASEP, the dynamics of the system will be determined 
entirely by the internal movement rates pi and p 2 . En- 
suring continuity of the current in a lattice with k p\ sites 
between each p2 site we find 



J = 



P\P2 



pI+V^Y 

P~2 



(T 2 



Pi + VP2 
Pi 



(2) 



Pi + y/P2 



These relations show the expected invariance under the 
interchange of p\ and p 2 . To find the boundary limited 
currents in the simple mean field approximation we as- 
sume that at least one of a or (3 is less than the smaller of 
Pi and p 2 . For the sake of simplicity, we also assume the 
lattice is composed of a whole number of periods, so that 
the last lattice site is a P2 site. Focusing initially on the 
entry-limited case and applying the current continuity 
conditions we find, for arbitrary k 



J = a(l-<Ti) = Pi<7i(1-<7 2 ) = . . . =P2(Tk+l (1-CTl ) (3) 

Yielding the following densities and currents for the k = 1 



J a = a(l-<Ti) = pi(Ti(l-(T2) = P2CT2(l-CTl) 
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Similarly in the exit limited case we have 
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J = P2CTfc+l(l - CTl) =Pl<7fc(l - CTfc+l) = . . 



(1) 



In writing Eqn Q we have assumed that the densities 
((Ti) in the lattice have the same period (T = (k + 1)) 
as the movement rates in the lattice [2^. Moreover, as 
the two-rate TASEP preserves the particle-hole symme- 
try found in the standard TASEP, we expect for a lattice 
of size N that cr.; = (1 — un-i+i) 24]. Coupling this con- 
dition to the assumption of periodicity in the densities, 
we find that within a single period, at = (1 — <Jk-i+2)- 
Thus we can immediately eliminate half the variables in 
Q). While straight-forward to solve, Q yields increas- 
ingly unwieldy expressions for the current and the den- 
sities as k increases. As a result we only show the result 
for the k = 1 case: 
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(6) 

To determine the transitions between the maximal cur- 
rent and entry limited regions we equate the maximal 
current solutions generated by to the expressions for 
J a and Jp in J2J and JSJ to find the boundaries be- 
tween the entry /exit limited and maximal current phases 
(a*, /?*). For the k = 1 case this yields 



(7) 
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Finally equating the the current expressions in equa- 
tions © and (JSJl we find that the transition between the 
entry and exit limited regions occurs when a = j3. 



B. Refined Mean Field Methods 

While appealingly transparent, we will see in the next 
section that the results of the simple mean field approach 
do not provide a particularly good match to the results 
of Monte Carlo simulations. 

We attribute the poor performance of the simple mean 
field method to the method's failure to capture the corre- 
lations between the occupation probabilities of different 
lattice sites. To address this deficiency we apply two al- 
ternative mean field approaches. The first refined mean 
field, which we call the finite-segment method (FSMFT), 
involves exactly solving the master equation for a fi- 
nite segment of the dual-rate TASEP lattice in a self- 
consistent manner (see Appendix . While primarily 
numerical in nature, this method can produce estimates 
for the current and density comparable to the results 
from MC simulations. The quality of the estimates pro- 
duced by the FSMFT approach is a function of the length 
of the finite-segment, with longer segments producing su- 
perior results. 

The second enhanced mean field method is an exten- 
sion to the two hopping rate TASEP of the well-known 
cluster mean field approach described in [t| and |15| . This 
method tries to capture the occupancy correlation be- 
tween adjacent sites by solving the exact master equa- 
tion for the multi-site probabilities P(ri, t^+i, ■ • ■)• 
As before Ti refers to the instantaneous occupancy of site 
i (see Appendix |B|) . As is typically the case, the cluster 
mean-field method becomes increasingly unwieldy as the 
size of the cluster increases. The results for the k = 1 
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As expected, the solutions are invariant under the ex- 
change (<7i,pi) <-* (02,^2), and regenerate the standard 
TASEP results in the limit p\ = p 2 . We note that in the 
maximal current solution, g\ + a 2 = 1 as a consequence 
of the particle-hole symmetry in the model. 

Finally, note that we can use these results to predict 
the location of the phase boundary between the maximal 
current and entry limited phases. Using the results for 
the densities in the maximal current phase and applying 
current continuity at the entrance of the lattice we find 
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To find the critical value of the extraction rate /?*, we 
enforce current continuity at the exit of the lattice 



/?* = 



2p\ + P2 
3(pi +P2) 

3ffllP2 

2(2 Pl +p 2 ) 



P1P2 



2(pi +P2) 
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The fact that a* = (3* is a reflection of the particle-hole 
symmetric intrinsic to the two-rate TASEP model. 

A similar approach (see Appendix 0) can used to find 
current and density approximations in the parameter 
regimes where the entry or exit rate limits the rate of 
particle transport through the lattice. First addressing 
the entry limited region we find 



a(pi(P2 ~ a) - c? + ^<jpl(p2 - a)a + (pip 2 - pia + a 2 ) 2 ) 
2p 2 (a+pi) 
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Similarly, the solution in the exit-limited phase is 
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(12) 



III. MONTE CARLO SIMULATIONS 

Extensive Monte Carlo simulation were performed to 
validate the various analytical models presented in the 
previous section. As we expected the densities in the 
lattice to vary significantly as the internal hopping rates 
were varied, we decided to base our Monte Carlo code 
on the BKL continuous-time algorithm 0|. The BKL 
algorithm has the advantage of maintaining a constant 
computational efficiency over a wide range of particle 
densities. 

The magnitude of the finite size effect in our sim- 
ulations was estimated by running lattices of varying 
lengths. For lattices at least one thousand sites long, 
the MC results were found to systematically deviate from 
the known TASEP results (pi = p 2 ) by less than half a 
percent. As a result, we used lattices containing approxi- 
mately 1000 sites for all our simulations. The simulations 
were run for 4 x 10 9 steps, which was sufficient to repro- 
duce the known TASEP results in lattices much longer 
than our 1000 site standard. In all simulations, p\ was 
normalized to 1 . Finally, to ensure an unbiased sampling 
of the lattice states, a linear-congruential pseudorandom 
number generator with a period of 2 x 10 was used [13. 



A. Currents 

The maximal current results (Figs [21 CD) show the ex- 
pected qualitative behavior with a low current for small 
values of p 2 and a value of 1/4 when p 2 = p\ = 1. In Fig 
|2]thc boundary-limited current predictions fEcin lll|) were 
used to generate a prediction for the maximal current as 
a function of p 2 . This was accomplished by maximizing 
J in Ecin llll with respect to a. Consistent with our ex- 
pectations, we find that the two cluster MF approaches, 
and the finite-segment mean field approach, provide bet- 
ter approximations to the Monte-Carlo results than does 
the simplest mean-field method. The relatively poor per- 
formance of the simple mean field model can be ascribed 
to the strong pair correlations present in the two-rate 
TASEP, which the simple MF model completely fails 
to capture. The 2-cluster maximal current mean field 
clearly performs best when p 2 /p\ ~ 1. In this, limit 
the model is essentially the normal single rate TASEP, 
and the expected pair correlations between adjacent p\ 
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FIG. 2: Current in the maximal current phase for a T — 2 
dual-rate TASEP. The results from the 2-cluster boundary 
limited MFT were produced by maximizing J in QUI with 
respect to a. 



and p 2 sites are negligible. As simple mean field meth- 
ods are known to provide exact results for the standard 
TASEP, the ability of the 2-cluster MF to provide accu- 
rate results in this region is not surprising. In contrast 
the boundary-limited 2-cluster MF performs best when 
P2 <C pi where the correlations between adjacent p\ and 
p 2 sites are expected to be strongest (although the known 
TASEP results are returned with p\ = p 2 ). Nonetheless 
in the limit p 2 — > 0, the relative error increases for all the 
mean field methods. 



B. Densities 

Referring to Figs [S] through \7\ we find the finite- 
segment MF, and both 2-cluster mean field methods pro- 
vide excellent matches to the simulation results. The 
densities produced by the maximal current 2-cluster 
MFT, the boundary limited 2-cluster MFT, and the finite 
segment method are all within 5% of the MC results, with 
the FSMFT results improving with increased segment 
length. For all three approaches the quality of the agree- 
ment is relatively uniform over all the values of p 2 , though 
the maximal current 2-cluster mean field results do tend 
to have an increased relative error at smaller values of 
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FIG. 3: (a) The effect of increasing segment length on the 
quality of the FSMFT estimates of the maximal current in 
a T = 2 dual-rate TASEP. The points are the Monte-Carlo 
simulation results, the lines are the FSMFT estimates. For 
segment lengths greater than or equal to twice the period, 
the quality of the current estimate shows only marginal im- 
provement with increasing length, (b) Maximal currents for 
dual-rate TASEPs of various periods. The FSMFT estimates 
were produced using a single period of the lattice as the finite 
segment. While providing reasonable estimates, the deviation 
between the MC and FSMFT results increases as the period 
(T) increases. As seen in (a), this increase in error can be 
mitigated by increasing the number of periods included in 
the finite-segment. 




FIG. 4: (a) Current profiles along the a direction in the T — 2 
dual-rate TASEP phase plane. The solid lines were produced 
using 11 U , the dashed lines were produced using (3J , and the 
dotted-dashed lines were produced using an N = 8 FSMFT 
approach, (b) Current profiles for dual-rate TASEPS of vari- 
ous periods with p 2 = 0.25. The FSMFT results (solid lines) 
were produced using an N = T FSMFT approach. From the 
plot we can observe the second order character of the entry- 
limited to maximal current transition. 



P2 where the correlations between adjacent sites are ex- 
pected to be strongest. In contrast, the simple mean field 
approach shows very significant deviations from the ex- 
perimental results particularly the small p2 limit. Given 
the simple mean field method's inability to capture cor- 
relations in the density between adjacent fast and slow 
sites, the poor performance of this method at small val- 
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FIG. 5: (a) Densities from the simple MF method, the 2- 
cluster mean field method, and the Monte Carlo (MC) simu- 
lations for a, T — 2 dual-rate TASEP. The simple mean field 
assumption results in the largest error for small values of p2 
where we expect the density correlation between adjacent sites 
to be large, (b) Results produced by the FSMFT method for 
various segment sizes. 
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FIG. 6: The bulk density profiles for T = 5 and T = 9 dual- 
rate TASEPS with P2/P1 = 0.2, in all three current phases. 
The plots show the good agreement between the Monte-Carlo 
(points) and N = T FSMFT (solid lines) results. 



C. Correlations 

Looking at the maximal-current density correlations 
displayed in Figs[S]and[51 we see that the pair correlation 
(71T2) is very large for small values of P2 and approaches 
zero as P2 — * Pi ■ The size of this pair correlation for small 
P2 confirms our intuition that the slow sites dominate 
the dynamics of the system, and is consistent with the 
maximal current assumption. 

Figure E{b) indicates that in the maximal current 
phase the (T2T1) pair correlation is negligible. The ab- 
sence of significant correlation in Fig|5Jb) and the rapid 
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FIG. 7: The density profile along the a direction in the T = 2 
dual-rate TASEP phase plane. The solid lines were produced 
using 1111 . the dashed lines were produced using |0J, while 
the dotted-dashed lines were produced by an N = 8 finite- 
segment MF approach. Note that although solution for 02 
is the same in Eqns. J1J and Hill , the predicted value of a* 
differs between the two mean field theories. This is the reason 
for the significant difference between the 02 profiles predicted 
by the two mean field methods at small values of p2 . 



FIG. 9: The pair correlation function produced by the MC 
simulations for the T — 2 dual-rate TASEP. Plot (a) dis- 
plays the correlation between fast sites and their neighbors 
to the right, while (b) does the same for slow sites. In both 
plots P2 = 0.1 In (a), the pair correlation (h — 1) is quite 
strong indicating that the slow sites are acting to limit the 
rate of transport through the lattice consistent with the max- 
imal current assumption. Conversely, (b) shows that slow 
sites are essentially uncorrelated with their neighbors to the 
right. 
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FIG. 8: The pair correlation function produced by the MC 
simulations and the various MF approximations for the T = 2 
dual-rate TASEP. As expected the 2-cluster mean field ap- 
proach provides a very good match to the MC results. As 
expected, the results produced by the segment MF improve 
with an increasing segment length. 
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FIG. 10: The phase diagram for the T = 2 two-rate TASEP. 
The dashed lines are produced by the simple mean field theory 
(Eqns 2] 03 ■ The dotted lines are produced by the maximal 
current 2-cluster method (Eqn|8J, and the solid lines are pro- 
duced by the boundary-limited mean field method (Eans tm 



decay of the correlation function in Fig Ela) indicates 
that the dynamics of the two-rate TASEP are essen- 
tially the dynamics of a set of successive, and nearly 
independent (pi,y>2) unit cells. While the correlation 
functions in Fig. [5] show longer range behavior, sites 
separated by more than two intervening lattice sites are 
essentially uncorrelated. Despite the maximal current 
2-cluster method's ability to provide highly accurate es- 
timates for the pair correlation function (Fig. |SJ, the 2- 
cluster method generates predictions for the current and 



density that carry small but significant error particularly 
for small p2- To seek the source of this error, we must 
consider the correlation profile produced by the maximal- 
current 2-cluster method. If we consider two lattice sites 
separated by h intervening lattice sites, the general pair 
correlation function can then be formed as follows: 
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r(h) = (TiT i+h ) / (n) (r i+1 ) - 1 

Where : 

, . _ Ex, [U^P(T it T i+1 )] Pfa^, 1) 

{nn+h) - u h - 1 P(n- -r \ 

L,{ Ti }^i=i P{n,n+i) 

(13) 

Referring to Fig. El we see that while the 2-cluster 
method provides good estimates for the pair correlation 
function, its predictions for the correlation profile are rel- 
atively poor. While the correlation profile produced from 
the results of the MC simulations decays rapidly to zero, 
the results are not negligible for h = 2,3, and thus we 
expect these small longer range correlations to have some 
effect on the current and density of the two site TASEP. 
The failure of the maximal current 2-cluster method to 
capture these weak longer range correlations may be re- 
sponsible for the increasing error in the mean field current 
and densities with decreasing p2 ■ 



D. Phase Diagram 

The phase diagram derived from simulations is dis- 
played in Fig along with the phase boundaries pre- 
dicted by the various mean field methods. The Monte 
Carlo phase boundaries were computed by taking numer- 
ical derivatives of the current as a function of a and (3 and 
locating any clear discontinuities in the derivatives. The 
dual-rate TASEP retains the general form of the standard 
TASEP phase-diagram, and the order of the transitions 
remains unchanged; first order between the high and low 
density regimes, and second order between the boundary 
limited and maximal current regions. The most signifi- 
cant deviation from the standard phase diagram is the in- 
crease in the maximal current region which accompanies 
a decrease in one of the hopping rates (j>i in our exam- 
ple) . Physically, the maximal current region is defined as 
the region of the (a, (3) parameter space where the inter- 
nal motions determine the net particle flow through the 
lattice. The increase in the area of the maximal current 
phase with decreasing p2 is then expected as we natu- 
rally require a or f3 to be less than p2 in order to enter a 
boundary limited region. 

Despite the varying degrees of success in predicting ac- 
curate steady state currents and densities, all three mean 
field approaches predict the phase boundaries within ap- 
proximately 12% error. As was observed in the mean field 
predictions for the steady state currents and densities, 
the boundary-limited 2-cluster mean field approach pro- 
vides the best predictions at small values of Pijpx-, while 
the maximal current mean field method excels when 
P2 ~ pi- We also note that the phase transitions pre- 
dicted by the boundary-limited 2-cluster approach were 
determined numerically by maximizing J in Eqns 111|) 
and H2J) 



IV. CONCLUSION 

We have developed and presented four separate ap- 
proximate methods to compute the the current and den- 
sities of totally asymmetric exclusions processes involv- 
ing two internal hopping rates. Moreover, we have ex- 
tensively simulated the two-rate TASEP and have thor- 
oughly explored its phase diagram. We find that the 
dual-rate TASEP retains the three phases found in the 
standard TASEP model. Within each of these phases 
the best of our mean field theories provide highly accu- 
rate, albeit not exact, predictions for the particle currents 
and densities. In particular, a maximal current phase 2- 
cluster mean field theory was developed that provides ac- 
curate estimates (Eqn. [SJ) for the current and the density 
(Figs. I2I5[) and the pair correlation function {t\T\) (Fig. 
|SJ. Similarly, in the boundary limited phases a simpli- 
fied 2-cluster MF method was introduced (Eqn. ^2 H2fl , 
which produced accurate estimates for the boundary cur- 
rents, and reasonably accurate estimates for the maximal 
currents (Fig. 0}. In addition to the generally inferior 
simple mean field approach (Eqns. IUEIII3), we have de- 
veloped a primarily numerical mean field method that we 
have termed the finite segment mean field method. This 
method can rapidly produce accurate estimates for the 
currents and densities in all three phases of the dual-rate 
TASEP. Additionally, the finite-segment method shows 
promise as a means of producing accurate estimates for 
the current and the density of TASEPs with hopping 
rates that vary spatially in a manner more complex than 
the periodic pattern described in this paper. 

Given the success of mean field theories in predicting 
the current and density in the standard TASEP, as well as 
various extended TASEPs, the inability of our mean field 
approaches to provide exact results is noteworthy. Com- 
parison of our best mean field predictions of the density 
correlation function with the Monte Carlo results (FigsOO 
[5J| indicate that while the 2-cluster mean field approach 
accurately captures the strong correlation between ad- 
jacent fast and slow lattice sites, the method does not 
account for the weaker, but still significant, longer range 
correlations. An extension to a three or four site cluster 
model may yield better analytic results. 
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APPENDIX A: THE FINITE-SEGMENT MEAN 
FIELD METHOD 

Here we give a detailed description of the the algo- 
rithm underlying the finite segment mean field method, 
and focus on using the method to solve for the currents 
and densities in the maximal current phase. Consider a 
finite lattice containing N sites and a whole number of 
(pi,Pi, ■ ■ ■ ,p%) periods. Working under the assumption 
that the bulk densities in the lattice are periodic and 
inherit the same periodicity as the movement rates, we 
focus on a single period. We then define the effective 
rate of injection into this finite segment of the lattice 
and the effective rate of extraction in terms of the un- 
known particle densities inside the segment. This yields 
"off = P2&k+i an d f3 c s = ^2(1 — 01 )■ Using these rates 
and the known values of p\ and P2 we can build the tran- 
sition matrix for the master equation describing the mo- 
tion of the particles in our T-site finite-segment. Letting 
N be the number of sites in the finite segment, and Ri 
be the movement rate for site i in the segment, the algo- 
rithm to construct the master equation transition matrix 
is (Fig. EJ: 





FIG. 11: The steps in the algorithm to generate the transition 
matrix for a TASEP model. For the purposes of illustration 
a three site model has been used. To begin we treat each 
possible occupancy of the lattice as a bit pattern, and label 
the state with the corresponding decimal value. For example 
given that Oil is the binary representation of 3, we label the 
state with particles occupying the second and third lattice 
sites as state 3. Next we divide the states into two groups; 
states where the first lattice site is occupied (1 -states), and 
states where the first lattice site is empty (0-states). Regard- 
less of the number of lattice sites in the TASEP, the transi- 
tions between the two classes of states always occur between 
the first half of the 1 -states and the second half of the 0-states 
(the straight green arrows in the figure). To determine the re- 
maining transitions we call the algorithm recursively on both 
the 1 -states and the 0-states, making sure to ignore the high- 
est order bit (i.e. the leftmost lattice site). Finally we add in 
the trivial transitions between each 0-state and each 1-state 
that are the result of injection at the left edge of the lattice 
(the red arrows in the diagram) . With the connectivity of the 
states fully enumerated, assigning the appropriate rate pi to 
each of the transitions is straightforward. 



1. Label every state of the TV-site finite segment with a 
number determined by the occupancy of the lattice 
sites. For example, an TV = 3 lattice with only the 
last site filled (ie an occupancy of 001), is said to 
be in state 1, while the same lattice with the last 
two sites filled (an occupancy of Oil) is said to be 
in state 3. 

2. Separate the states generated in step ^ into two 
groups; the states where the left-most lattice site 
is empty (the 0-states), and the state where the 
left-most lattice site is full (the 1 -states). 

3. Enumerate the transitions between the 1-states and 
the 0-states. The transitions will always be the 
movements that take the system from the lowest 
2 N ~ 2 1-states to the highest 2 N ~ 2 0-states. The 
rate of these transitions is Ri , where I is the current 
depth of the recursion (see stepBJ. 

4. Recurse on the 0-states and then on the 1-states, 
ignoring the leftmost site (ie with each recursion of 
the algorithm, we reduce the length of the segment 
by 1), or stop if the current segment length is N = 
1. 

5. Once the recursion specified in steps 2] through 0] 
has completed, enumerate all the trivial transitions 
between the 0-states and the 1-states that occur 
at rate a (ie the state transitions due to particle 
injection on the left). 

With the master equation constructed by this algo- 
rithm, we can then compute the densities in the finite 
segment, apply a self consistency condition on the den- 
sities, and subsequently compute the current in terms of 
the hopping rates p\ and p2- Thus utilizing the state 
enumeration scheme and algorithm described in Fig. 1111 
we generate and solve the following equations 



M({a i },p 1 ,p 2 )V t 



= 



-iV i ({a i },p 1 ,p 2 ) = 01 



(Al) 



Ei=l V 2 i-l({cri},p 1 ,p 2 ) = (T N 



In Eqn IjAljl . M is the transition matrix, while V{ is 
the probability of finding the finite-segment in occupancy 
state < i < (2 N - 1) (see Fig. El)- When applying the 
finite-segment mean field approach, I jAlj l is solved itera- 
tively until a fixed point in the densities is reached. Using 
the results from IjAljl we can compute the current from 

2N-1 2 JV_1 

the expression J = a c g J2i=i v i ( or &ffI3i=i v 2i-i)- 
Additionally, the finite-segment method can be extended 
to treat the boundary limited dual-rate TASEP by fix- 
ing a e ff = a or /3 e gp = (3 as appropriate, and applying 
the density self-consistency condition at the end of the 
finite-segment that lies in the interior of the TASEP lat- 
tice. While a simple and effective approach, the exponen- 
tial increase in the size of the transition matrix with the 
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length of the segment quickly renders IjAlfl analytically 
untreatable. As a result we chose to solve (IA1|) numer- 
ically using a simple fixed-point iteration scheme. Em- 
ploying the fast Arnoldi-method eigensolvers in the well- 
known ARPACK linear-algebra software library 0] , we 
could easily treat segments containing ~ 20 lattice sites. 



dP(Tj,T i+1 ) _ 

dt ~ 

-Q(Ti-l,Ti,T i+1 ,T i+ 2)(pi5 Titl S Ti+u o) 

+ Et,_ 1 t,_ 2 P(n-2,n,T i - 1 ,T i+1 )(p 2 S T ^ 1 5 Ti _ u0 ) 



APPENDIX B: CLUSTER MEAN FIELD +L«r i , P(n,T. l+2 ,T l+1 ,T l+3 )(p 2 S taUl+u0 5 Tt+2yl ) 



METHODS 



~P(Ti, T i+ i, T l+2 , T i+3 )(piS n+1 ^S Ti+2 fi) 



We begin the derivation of the current and density ap- 
proximations displayed in Eqn (|SJ by considering the two 
site probability P(ri,T i+ i). Here, refers to the occu- 
pancy of site i, and is equal to 1 if the site is occupied, 
if the site is empty. The time evolution of the occupancy 
state of any two adjacent sites in a TASEP will depend 
on the two sites themselves along with the site immedi- 
ately to the left and the right of the two site grouping 
(FigEb)- 



(Bl) 

As we have incorporated two different hopping rates 
into the TASEP model we are forced to deal with two sets 
of probabilities; the first (P(rj, Ti+i, Tj_|_2> Tj+3)) refers to 
the probability of having configuration (t^, Tj-|_i, . . .) given 
that the first site (site i) is a p\ site. The second prob- 
ability term (Q(tj, Tj+i, Ti+a> Ti+3)), is the probability of 
configuration (r,-, Tj+i, . . .) when the first site is a p 2 site. 
Now we assume we can factor the 4-site probabilities as 
follows 



P 
'2 



P 
'2 



Tft 



eff 



N = 4 Segment 



P 2 



P 2 



R 



R 



R 



+ 1 



FIG. 12: The lattice site groupings used in the two enhanced 
mean field approaches, (a) The arrangement used for an N = 
4 finite segment mean field (MFT) approach. The master 
equation for the state of the four marked lattice sites is solved 
exactly. The four sites are coupled to the rest of the lattice 
by assuming an effective injection rate of a e s = on the 
left, and an effective extraction rate of /? e ff = P2(l — 01) on 
the right, (b) The arrangement of lattice sites used in the 
2-cluster mean field theory. The two lattice sites of interest 
are the central two sites in grouping Rq. Groupings J?+i and 
R-\ are produced by shifting the boundaries of the group 
Ro one lattice site to the right, and one lattice site to the 
left respectively. To generate a second set of equations, we 
simply interchange all the pi and P2 sites, while keeping the 
groupings R+i, Ro, and the same. 



Thus the master equation for the two site probability is 



P\n,T i+ l,T l+2 ,T i+3 ) = P(Ti,T l+ l)Q(T l+ l,Ti +2 )P(T l+2 ,Ti +3 ) 
Q{l~i,T t +2,T l+2 ,Ti +3 ) = Q(T l ,rj + l)P(rj + l,Tj + 2)Q(rj +2 ,T l+3 ) 

(B2) 

And by definition we have 



o-i = P(1,0) + P(1,1) 
a 2 = Q(1,0) + Q(1,1) 
P(0,0) = 1-P(0,1)-P(1,0)-P(1,1) 

Thus, in the steady state we can write 
dP(0,0) 



dt 



P2 (P(0,1,0,1)-P(1, 1,0,0)) 
p 2 (P(0,l)Q(l,0)P(0,l) 

-P(1,1)Q(1,0)P(0,0)) = 



Upon using i|B3|) . we find 

0-i(l - <J 2 ) - (cr 2 - <7i) 2 



P(1,0) 



1 + ((72 - CTi) 



(B3) 



(B4) 



(B5) 



The equation for Q(l, 0) can be found from l|B5|) . by in- 
terchanging pi and p 2 , as well as ai and a 2 . Applying 
the current continuity condition piP(l,0) = p2Q(l,0), 
we then arrive at the relation 



Pi- 



Qi(l - 02) - {02 - gi) 
1 + (cr 2 — (Ti) 



P2- 



cr 2 (l - CTi) - (01 - <7 2 ) 2 



1 + (ai - a 2 ) 



(B6) 



10 



We now make the substitution a% — 
Pi to be the larger of the two rates. 



ax = fc and define 
Thus Hi < er 2 and 



J 



fc > 0. Solving i|B6(l for ai, we find o\ 



[pa(l - fc 2 ) - pi(fe - l) 2 ] + ^(1 - fc 2 )(fc(pi + p 2 ) - (pi - P2))(3A;(pi + p 2 ) - (pi - p 2 )) 



2/c(pi + p 2 ) - (pi 



r 



P2j 



(B7) 



From (|B7|I we find that <j\ is real only when k < 

or fc > ■ Substituting (JEJ into JElJl we find J = 

7 k ^lf , . , which shows J < for fc > 21^2= J > 

for < fc < and U = piP2(pi-P2) > Qyk 

P1+P2' rffe (fe(P2+Pl) + (P2-pl)) 2 

Since J must be positive, ci must be real, and J is mono- 
tonically increasing, the maximum value of the current 
must occur when fc = „f 1 ~ P2 \ . Substituting this value 

3(Pl+P2) b 

of fc into (|B7p . we arrive at ijBJl. 

An analogous approach can be used to find current and 
density estimates in the entry and exit limited regions. 
Considering the entry limited case first, we find the fol- 
lowing master equation for the occupancy of the first two 
sites in the lattice: 



In writing (|B8|I we have assumed that the occupancy 
probabilities P(t\, r 2 ) do not vary significantly as a func- 
tion of position near the entrance of the lattice. This 
assumption is analogous to that used with the standard 
TASEP in the boundary limited region where the den- 
sity is assumed to be essentially constant. Applying the 
normalization condition on the probabilities P(Ti, Tj_|_i) 
and solving (|B8p generates the results displayed in lUfl . 
An analogous approach yields the current and densities 
in the exit limited regime (Eqn. I12|) . 



dP(0, 0) 
~dt 

dP(0, 1) 
dt 



dP {1,0) 
dt 



dt 



-aP{0, 0) + P2 P(0, 1)(P(0, 0) + P(0, 1)) 

-aP(0, 1) - P2 P(0, 1)(P(0, 0) + P(0, 1)) 
+PiP(l,0) 

aP(0, 0) +p 2 P(l, 1)(P(0, 0) + P(0, 1)) 
-PiP(l,0) 

- P2 P(1, 1)(P(0, 0) + P(0, 1)) + aP(0, 1) 

(B8) 
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